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We study the one-band Hubbard model on the trellis lattice, a two-dimensional frustrated lattice 
of coupled two-leg ladders, with hopping amplitude t within ladders and t' between ladders. For 
large U/t this is a model for the cuprate Sri4_a;Caa;Cu2404i. We investigate the phase diagram 
as a function of doping for U = lOt using two quantum cluster methods: The variational cluster 
approximation (VGA), with clusters of sizes 8 and 12, and Cellular dynamical mean field theory 
(CDMFT), both at zero temperature. Both methods predict a superconducting dome, ending at 
roughly 20% doping in VGA and 15% in CDMFT. In VGA, the superconducting order parameter 
is complex in a range of doping centered around 10%, corresponding to bulk chiral, T-violating 
superconductivity. However, the CDMFT solution is not chiral. We find evidence for a migration 
of the Cooper pairs from the inter-ladder region towards the plaquettes as doping is increased. 


I. INTRODUCTION 

Inspired by the discovery of high-Tc superconductivity, 
Dagotto et al. predicted a superconducting phase in 
a theoretical model consisting of weakly-coupled, quasi- 
one-dimensional ladders. This model exhibits a spin-gap 
and d-wave hole-pair formation away from half filling [1- 
3] . This prediction was realized in the copper oxide-based 
ladder material Sri4_a;Caa;Cu2404i with hole doping. At 
X = 13.6, the critical temperatures under pressures of 3 
GPa and 4.5 GPa are Tc = 12K and Tc = 9K, respec¬ 
tively [4]. At x = 11.5, a superconducting dome is seen as 
a function of pressure [5]. Recently, the critical tempera¬ 
ture of the a; = 11 compound Sr 3 CaiiCu 2404 i has been 
doubled, from 12K to 24K, by applying a weak uniaxial 
pressure of 0.06 GPa [6]. 

Many theoretical studies have been reported on the 
single-ladder Hubbard model, with and without doping, 
using a variety of methods: Exact diagonalizations [2, 7], 
density-matrix renormalization group [8], resonating- 
valence-bond (RVB) mean-field theory [9], bosoniza- 
tion [10] and quantum Monte Carlo (QMC) [11, 12]. A 
consistent result from those studies is the emergence of d- 
wave superconducting correlations in the double leg lad¬ 
der upon doping. Coupled ladders described by the trel¬ 
lis lattice have been investigated using the Fluctuation 
Exchange (FLEX) method, confirming the possibility of 
d-wave superconductivity at half-filling [13]. 

In this paper, we report on a theoretical study of 
superconductivity in the one-band Hubbard model on 
the trellis lattice away from half-filling, at zero temper¬ 
ature. We use the Variational Cluster Approximation 
(VGA) [14] and Cellular Dynamical Mean-Field Theory 
(CDMFT) [15, 16]. For the range of on-site repulsion 
studied, superconductivity does not occur at half-filling, 
but a superconducting dome appears upon doping. More¬ 
over, the superconducting order parameter computed 
from VGA becomes complex in a range of doping centered 
around 10%, thus breaking time-reversal symmetry. The 
energy gain from the chiral nature of superconductivity 


is small, at best of the condensation energy, and the 
chiral solution is not found with CDMFT. 

This paper is organized as follows. In Section H the 
model is presented, as well as the structure of singlet 
superconductivity in the BCS approximation. In Sec¬ 
tion HI the VGA technique is summarized and the results 
of its application are presented; this is the main part of 
the paper. In Section IV CDMFT is applied in order to 
confirm by an independent method the occurrence of su¬ 
perconductivity. We add a short discussion and conclude 
in Section V. 


FIG. 1. (Color online) Model for the ladder cuprate 
Sri4_a;Caa;Cu2404i ou the trellis lattice. The shaded areas 
are the clusters used in VGA, labeled C8 and C12. 


II. MODEL AND MEAN-FIELD 
REPRESENTATION 

The Hamiltonian of the one-band, repulsive Hubbard 
model on the trellis lattice is 
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where Cia annihilates an electron of spin a at site ty = 
tji is the hopping amplitude between sites i and j, riia = 
cJ^Cicr is the number operator at site i and U the on-site 
Coulomb repulsion. The density of electrons is controlled 
by the chemical potential fi. The only nonzero hopping 
terms are indicated by inter-site links on Fig. 1, with 
tij = t on the ladder plaquettes, and Uj = t' between 
ladders. Model (1) has two bands: the unit cell contains 
two orbitals, separated by e^, on the figure. The vectors 
eo and on Fig. 1 define a basis for the lattice. For 
convenience, we will define the Brillouin zone exactly like 
on the graphene lattice. 

This model offers an approximate description of 
Sr 3 CaiiCu 2404 i, wherein each site represents a copper 
atom. In the actual material, oxygen atoms are located 
midway between copper atoms on the links of each square 
plaquette. The relation between hole doping 5 (the elec¬ 
tron density is n = 1 —(5) and Ca doping x in the material 
is not simple, as 5 is also affected by pressure. Through¬ 
out this paper we will set t = 1 and t' = 0.15; this ratio 
t'/t is taken from band structure calculations [17]. The 
value oiU/t will be set to 10 in most VGA and CDMFT 


/-2cos(k • ^ 7k 

7^ -2cos(k-ej^) 

V < ^k 


with 7 k = —. This is the most 
general form for singlet superconductivity. If we assume 
only nearest-neighbor pairing with amplitudes Dx^ —Dy 
and Di 2 in the directions e^,. By and ei _2 respectively, 
the anomalous terms of that matrix are 

0k = —Dy COS(k • By) 

ryk = + D2e**^"== (5) 

The choice of sign for Dy reflects our anticipation of d- 
wave superconductivity on the plaquettes. 

If the superconductor is time-reversal (T) invariant, 
the components of i7k satisfy the relation i7_k = 

This implies that the amplitudes 7?a;,y,i,2 defined above 
are all real. On the other hand, if any one of them is 
complex, the superconductor breaks time-reversal invari¬ 
ance. 

Let us stress that we are not performing a true mean- 
field computation here: there is no factorization of the 
interaction, no self-consistent procedure, etc. Indeed, the 
Hubbard model, with its local repulsion, is not amenable 
to a self-consistent (BCS) mean-field computation of d- 
wave superconductivity. We are simply illustrating the 
form that superconductivity can take in this model, in 
order to compare with the complete variational or self- 
consistent computations reported on in Section III. 


computations. 


A. Superconductivity 

What form can superconductivity take in such a 
model? To answer this question, let us first provide a de¬ 
scription of the superconducting order-parameter at the 
mean-field level. It is then convenient to adopt a Nambu 
description, with the multiplet of destruction/creation 
operators 

Ck = (ci^(k),c2t(k),cj^(-k),c^^(-k)) (2) 

where the first index of each operator is a sublattice in¬ 
dex, distinguishing the left and right sites of each rung. 
This array of operators is used in a matrix description of 
the non-interacting, mean-field Hamiltonian 

Hbcs = E 4i7kCk (3) 

k 

with the 4x4 Hermitian matrix 


0k Vk 

- M ??-k 0k 

2 cos(k ■ By) + fi - 7 k 

-7k 2cos(k • Gy)-I- 



B. Order parameter 

The most general way to represent superconducting 
order is to plot the momentum-dependent order param¬ 
eter Aab(k), defined as the integral over frequency of the 
Gorkov function (the anomalous part of the Green func¬ 
tion): 


/ OO J 

^Fo&(k,fw) (6) 

Here (a, b) are sublattice indices taking two possible val¬ 
ues, associated with the left and right legs of the ladder. 
The Gorkov function Fab is the top-right block of the 
Nambu Green function G^y(k, w) defined as follows at 
zero temperature: 


= (H|C^(k)-J_^ct(k)|H) 

LU — JGL -r h/Q 

+ (H|Ci(k) ^ C^(k)|H) , (7) 

where w is a complex-valued frequency, |H) is the many- 
body ground state and Eq the ground state energy. For a 
two-band model, Fab = Ga.b+ 2 - In the special case of the 
non-interacting BGS Hamiltonian (3), the Nambu Green 
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FIG. 2. (Color online) Superconducting order parameter 
Aa[,(k) for the BCS Hamiltonian (3), as a function of wave- 
vector, with parameters ~ 0.09, Dy = 0.03 and D\ = 
D 2 = —O.Ol-fO.OOSi. Doping is set at 10%. The Brillouin zone 
is indicated. Panels (a) and (b) show the real and imaginary 
parts of the rung component (different sublattices), whereas 
panels (c) and (d) show the leg component (same sublattice). 
Red means negative, blue positive and the range is indicated 
on top of each panel. The normal state Fermi surface is shown 
in red on panel (a). 


function is 

<«) 

In order to connect with the more familiar one-band 
BCS theory, let us point out that in that case the matrix 
iJk has the simpler form 




Ck - It ^'k A 
^k —Ek + At/ 


( 9 ) 


where Ck is the dispersion relation and 0 k the gap func¬ 
tion. The order parameter is then simply 


Ak 


0 k 

2\/ ^k + l^kp 


( 10 ) 


Figure 2 illustrates the superconducting order param¬ 
eter Aab(k) in a particular case of Model (3). This will 
later be compared to a solution, obtained through VGA, 
that contains correlation effects. We have chosen super¬ 
conducting amplitudes Dy, Di,D 2 ) that break time 
reversal slightly and match local order parameters of an 
actual VGA solution described later on. The fact that 
ImAi 2 (k) 7 ^ — ImAi 2 (—k) and ImAii(k) 7 ^ 0 is a vi¬ 
sual signature of time-reversal breaking. 


III. THE VARIATIONAL CLUSTER 
APPROXIMATION 

We use the Variational Cluster Approximation 
(VGA) [14] to investigate the zero-temperature phase di¬ 
agram of Model (1), more specifically the existence of 
d-wave superconductivity upon hole doping for several 
values of U. VGA - also called VCPT in its early days 
- has been used to study the emergence of d-wave su¬ 
perconductivity in a simple description of the high-T,. 
cuprates based on the Hubbard model [18, 19]. It is based 
on Potthoff’s self-energy functional approach [20]. For a 
review, see Ref. [21]. 


A. Description of the method 

In VGA, we must distinguish between the original 
Hamiltonian H, defined on the infinite lattice, and a 
reference Hamiltonian H', defined on a small cluster of 
atoms. H' is a restriction of H to the cluster, except 
that a finite number of Weiss fields may be added to it, 
in order to probe certain broken symmetries. Any one- 
body term can also be added to H'. The electron self¬ 
energy S(w) associated with H' is used as a variational 
self-energy, in order to construct the Potthoff self-energy 
functional: 

n[E{h)] = H'[S(/i)] 

+ Tr ln[-(Go 1 - S(d))-i] - Tr ln(-G'(d)) (11) 

where G' is the physical Green function of the cluster. 
Go is the non interacting Green function of the origi¬ 
nal model and h denotes collectively the coefficients of 
all the adjustable one-body terms added to H' acting 
as variational parameters. The symbol Tr stands for a 
functional trace, i.e., a sum over all degrees of freedom 
(e.g. momenta or sites) and frequencies. is the ground 
state energy (chemical potential included) of the cluster 
which, along with the associated Green function G', is 
computed numerically, in our case via the exact diago- 
nalization method at zero temperature. 

Eq. (11) provides us with an exact, non-perturbative 
value of the Potthoff functional r2[S(/i)], albeit on a re¬ 
stricted space of self-energies E(/i) which are the physical 
self-energies of the reference Hamiltonian H'. Expression 
( 11 ) is computed numerically in order to look for station¬ 
ary points of that functional, for instance via a Newton 
or quasi-Newton method. The resulting value of h de¬ 
fines the best possible self-energy S(a;) for that parame¬ 
ter set; the latter is then combined with Go to form an 
approximate Green function G for the original Hamilto¬ 
nian H, from which any one-body quantity, for instance 
the order-parameters associated with broken symmetries, 
can be computed. 

When confronted with competing solutions, i.e., differ¬ 
ent stationary points of r 2 [S(/i)] or points obtained via 
different sets of Weiss fields, the one with the lowest value 
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FIG. 3. (Color online) Pairing order parameters (A^, A^, A^) 
as a function of doping S for the two clusters used in VGA, at 
U = 10. The complex phase of A^ can be read on the right 
axis. The T-breaking solution (arg Ax 7^ 0) exists in a finite 
range of doping. 


of the Potthoff functional is selected, as is an approx¬ 
imation of the exact free energy of the infinite system. 
VGA retains the correlated character of the model, since 
the local interaction is not factorized. The approxima¬ 
tion may be controlled in principle by varying the size 
of the cluster and the number of variational parameters 
used. 

In this work we use the clusters labeled C8 and C12 
illustrated on Fig. 1. It is important to test more than 
one cluster, as there will be some variance in numerical 
results among different clusters and robust characteristics 
need to be identified. Larger clusters will generally lead 
to smaller values of the order parameter, because of in¬ 
cluster order parameter fluctuations. 


B. Superconductivity 


In VGA the possible presence of d-wave superconduc¬ 
tivity is probed by adding to the cluster Hamiltonian H' 
pairing operators. These may be expressed in terms of 
the singlet pairing operators — Ci^Cj^ = Kji 


as follows: 


Ax = Aij 

ri-rj=e,c 

A = - E 

ri-rj=ey 

Ai= E 

ri-rj=ei 

A 2 = \ 

ri-rj=e2 


where the unit vectors ex^y^i ^2 are defined on Fig. 1. A^, 
is the sum of all pairing operators along rungs. Ay is the 


sum of all pairing operators along legs, and Ai 2 are the 
sum of pairing operators between the ladders. In prac¬ 
tice, the Nambu description (2) is used: A particle-hole 
transformation is applied to spin-down operators, giving 
the above pairing operators the appearance of hopping 
terms. 

The Weiss Hamiltonian added to the cluster Hamilto¬ 
nian H' takes the form 


idgc — Dx^x + DyAy -b DiAi + D 2 A 2 + H.c. (13) 

where the coefficients {Dx, Dy, Di, D 2 ), the so-called 
Weiss fields, are variational parameters, adjusted so as 
to make the Potthoff functional (II) stationary (in prac¬ 
tice, minimum). The minus sign in front of Ay in (12) 
means that we anticipate d-wave symmetry on the pla- 
quette, i.e., we anticipate Dx and Dy to have the same 
sign, which is indeed what we find numerically. 

In principle, the Weiss fields may be complex-valued, 
as id'j. remains Hermitian anyway. However, we choose 
Dx to be real so as to fix the overall phase. Then Dy, Di 
and Z ?2 can be complex. In practice, in order to limit the 
number of variational parameters, we assume that Dy is 
real and that Re I ?2 = ± Re Di and Im D 2 = ± Im Di. 
We found that the lowest minima of the Potthoff func¬ 
tional have Di = D 2 = D±, and in the rest of this paper 
we will accordingly define Ax = Ai -|- A 2 , for a total of 
4 variational parameters: Dx, Dy, ReUx and Iml?x- 

Figure 3 shows the order parameters A^ = (Aq.) -b 
(Alj), for a = x,y,A, computed from the VGA Green 
function, as a function of hole doping, for C/ = 10 and 
the two clusters shown on Fig. 1. The inter-ladder or¬ 
der parameter Ax is complex over a range of doping: 
its modulus is plotted, along with its phase (right ver¬ 
tical axis). The superconducting dome has a maximum 
between 10% and 15% (depending on the cluster) and 
ends at about 20% doping. It falls to zero exactly at 
half-filling. The inter-ladder order parameter Ax is no¬ 
ticeably smaller than the ladder order parameters A^, and 
Ay, but this is roughly in line with the ratio t'/t = 0.15. 
The rung and leg order parameters (A^, and Ay) also 
have imaginary parts whenever Ax has one, but they 
are small and would not make visible contributions to 
|Aa;| or I Ay I on the plots. 

The order parameters shown on Fig. 3 are special con¬ 
volutions of the general momentum-dependent order pa¬ 
rameter Aaf,(k) with particular form factors associated 
with nearest-neighbor pairing. They have the advantage 
of simplicity, but are somewhat arbitrary. Unfortunately, 
the full order parameter Aa&(k) cannot be plotted sim¬ 
ply as a function of doping. However, Figure 4 shows 
Aa&(k) for the VGA solution at 10% doping. This is to 
be compared with Fig. 2, which shows the corresponding 
BGS order parameter, obtained by setting the BGS fields 
to values that reproduce the same values of the link or¬ 
der parameters Ax^y^x- We notice that the features of 
Fig. 4 are qualitatively the same as those of Fig. 2, al¬ 
though less sharp, because of strong correlation effects. 
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FIG. 4. (Color online) Superconducting order parameter 
Aa[,(k) computed from the VGA solution, as a function of 
wave-vector. 17 = 10 and doping is 10%. This is to be com¬ 
pared with the BCS order parameter of Fig. 2. 


C8,U= 10 



FIG. 5. (Color online) Loop supercurrents Ip and Ic circu¬ 
lating respectively around a plaquette and around the center 
parallelogram of the C8 cluster, as a function of doping, for 
U — 10. The loops are illustrated on Fig. 1 


The sharp lines of Fig. 2 have become broad maxima and 
minima, but the asymmetry of Im A 21 stands out. Note 
that the scales (color range) differ from those of Fig. 2 
by factors of two to three. 

The T-breaking nature of the solutions found can also 
be assessed by computing chiral supercurrents. Fig. 5 
shows the supercurrents Ip and Ic circulating along the 
loops indicated on Fig. 1. These are defined as the ex¬ 
pectation values of 

/ = V {4,<TCi+i,a - H.c.^ (14) 

2Gloop,(7 



FIG. 6. (Color online) Value of the Potthoff functional at 
the solution as a function of chemical potential /r for U — 10 
and different sets of variational parameters: the normal solu¬ 
tion (no variational parameters), two T-preserving solutions, 
with {Dx,Dy) and {Dx,Dy, ReDx) used as variational pa¬ 
rameters, and the T-violating solution with variational set 
{Dx,Ry, ReDx, ImUx). The latter has the lowest value of 
fl. Two values of doping are indicated for that solution. A 
multiple of /i is added in order to better separate the different 
curves. 


where the sum is taken around the loop. The expecta¬ 
tion values (/p,c)cius in the cluster ground state vanishes 
if Ax is real, but is nonzero as soon as Ax develops an 
imaginary part. This demonstrates that, in the latter 
case, superconductivity is chiral: if ImAx changes sign, 
the value of the Potthoff functional does not change - 
hence we again have a VGA solution - but the sign of 
the current changes. Note that these supercurrents are 
measured on the cluster itself, as ground state expecta¬ 
tion values, without using the Green function, because 
the latter provides expectation values on the whole lat¬ 
tice and these current loops cancel each other when the 
loops are stacked on the lattice. They are computed only 
to underline the chiral character of the complex super¬ 
conducting solutions. Note that the supercurrent loop 
located between the ladders dominates at small doping, 
whereas the contrary is true of the plaquette supercur¬ 
rent. This leads us to believe that the Gooper pairs tend 
to locate between the ladders at small doping and move 
towards the plaquettes at larger doping. 

The T-breaking solution has the lowest energy in a 
sequence of solutions that can be obtained in VGA by 
increasing the number of variational parameters, as il¬ 
lustrated on Fig. 6. We plot the value of the Potthoff 
functional at the solution, which is an approximation to 
the grand potential 12, as a function of chemical poten¬ 
tial /i, since 12 is by construction as a function of /r, not 
density. In the top curve, no variational parameters were 
used, and this represents the normal solution. The sec¬ 
ond curve from the top is obtained by using the plaquette 
anomalous Weiss fields {Dx,Dy) as variational parame¬ 
ters. The third curve is obtained by adding the real part 
of the inter-ladder pairing D± to the set. Finally, the 
lowest curve is obtained by adding both the real and 
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imaginary parts of D±^ to the set, and the correspond¬ 
ing solutions break time-reversal invariance. This illus¬ 
trates the process by which the quality of VGA solutions 
is improved by adding variational parameters. Another 
solution, obtained by allowing Dy to take complex val¬ 
ues, is not shown, as it is hardly distinguishable from 
the last one. We see that the energy advantage of the 
T-breaking solution is at best 0.0025t, or roughly A Qf 
the condensation energy (the difference between the nor¬ 
mal state and superconducting state energies), and this 
only at the most favorable doping {5 ^ 12%). Thus, even 
though the VGA simulation predicts T-breaking super¬ 
conductivity in this system in a range of doping, it must 
be kept in mind that this solution is very close in en¬ 
ergy to other approximate solutions that do not break 
time-reversal symmetry. 

Notice that the difference between the second and third 
curves on Fig. 6 increases with /r, i.e., towards smaller 
doping. This means that the importance of varying 
is greater on the underdoped side of the dome, which 
confirms our interpretation that the Gooper pairs tend 
to locate between the ladders at small doping. 

IV. THE CELLULAR DYNAMICAL MEAN 
EIELD THEORY 

We also used Cellular Dynamical Mean Field Theory 
(CDMFT) to confirm the appearance of a superconduct¬ 
ing dome by independent means. 

A. Description of the method 

CDMFT like VGA, proceeds by tiling the lattice with 
clusters and by computing an optimized self-energy for 
each cluster. Unlike VGA, the space of self-energies is 
not explored by adding Weiss fields on the cluster, but 
rather by coupling each cluster to a bath of uncorrelated, 
auxiliary orbitals that represent the effect of the cluster’s 
environment [15, 16, 22, 23]. The cluster Hamiltonian 
is supplemented by bath-cluster hybridization and bath 
energy terms: 

.^bath — ^ ^ -f ^ F H.C. (1^) 

where denotes the annihilation operator for the bath 
orbital labeled 

This, together with the restriction of the Hubbard 
Hamiltonian (1) to the cluster, defines an Anderson im¬ 
purity model. The cluster Green function, when traced 
over the bath orbitals, takes the following form as a func¬ 
tion of complex frequency ui: 

G'-^{uj)=uj-t-T{uj)-'E{uj) (16) 

where the hybridization matrix r(a;) is 

r(a;) = e(a; - e)-i6>^ (17) 



FIG. 7. (Color online) CDMFT cluster-bath system used in 
this work. The trellis lattice was deformed into a square lat¬ 
tice for simplicity. Right panel: arrangement of the two clus¬ 
ters A and B needed to tile the lattice. Left panel: picto¬ 
rial representation of the two clusters: lattice sites are gray 
spheres and bath orbitals are represented by smaller, red 
spheres. Green links represent hopping terms, dashed blue 
lines are bath-cluster hybridization terms and red dotted lines 
anomalous terms between bath sites, forming the anomalous 
part of the matrix e. 

in terms of the matrices and e^i,. The Green function 
G(k, w) for the lattice model is then computed from the 
cluster’s self-energy as 

G-i(k,cc) = Go-'(k,u;)-S(u;) (18) 

Here k denotes a reduced wave-vector, belonging to the 
Brillouin zone associated with the superlattice of clusters 
that defines the tiling. All Green function-related quan¬ 
tities are 2A(. x 2Nc matrices, W being the number of 
sites in the unit cell of the superlattice, which is made of 
one or more distinct clusters (the factor of 2 is there be¬ 
cause of spin). Go is the non-interacting Green function. 
In practice, the cluster Green function is computed from 
an exact diagonalization technique using variants of the 
Lanczos method (just like in VGA). Then the self-energy 
is extracted from Eq. (16). 

The bath and hybridization parameters {tyv,Qa^i) are 
determined by the self-consistency condition 

= (19) 

ic 

{N is the [quasi-infinite] number of sites in the whole 
system). In other words, the local Green function 
G'( Lo) should coincide with the zero wave-vector Fourier 
transform of the full Green function. This condition 
should hold at all frequencies, which is impossible in a 
zero-temperature implementation of GDMFT because of 
the finite number of bath parameters at our disposal. 
Therefore, condition (19) is only approximately satisfied. 
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FIG. 8. (Color online) Order parameters computed from the 
CDMFT solutions found with the cluster-bath system illus¬ 
trated in Fig. 7, as a function of doping. The CDMFT solu¬ 
tions do not break time-reversal invariance. 



0 0.05 0.1 0.15 0.2 
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through the use of a merit function. Details can be found, 
for instance, in Ref. [23]. 


FIG. 9. (Color online) Order parameters (top) and |Aj_| 
(bottom) as a function of doping 5 for cluster C8 and several 
values of on-site repulsion U. 


B. System studied and superconductivity 

When modeling superconductivity in CDMFT, it is 
convenient to introduce anomalous terms between bath 
sites, thus treating bath sites as if they were forming a 
‘phantom cluster’. The cluster-bath system used in this 
work is illustrated on Fig. 7. Two unequivalent, four-site 
clusters form the repeated unit cell of the super-lattice. 
The Nambu formalism is used to represent anomalous 
terms (see Ref. [24] for explanations of its use in the con¬ 
text of CDMFT). The bath orbitals are grouped into two 
sets of four, and within each set anomalous terms are de¬ 
fined that mimic what could occur on the cluster itself 
(hence the expression ‘phantom cluster’). Each bath set 
has four or five links (dotted red lines on the figure) and a 
complex pairing operator is defined on each of these links, 
except on the rung link where it is assumed to be real, in 
order to set the global phase of the superconducting state. 
Taking symmetries into account, this makes for a total 
of 14 bath parameters for superconductivity, in addition 
to 8 bath orbital energies and as many hybridization pa¬ 
rameters, for a total of 30 variational parameters. 

Figure 8 shows the order parameters A^,, Ay and A± 
for the CDMFT solutions obtained at [/ = 10. These so¬ 
lutions do not break time-reversal in any significant way 
(Im Ax < 10“®). But the superconducting dome seen in 
VC A is still there, although somewhat narrower (noth¬ 
ing beyond <5 = 15%). As doping is increased, the inter¬ 
ladder pairing operator Ax has a first maximum around 
1.5%, then decreases before increasing again, carried by 
the other components. This is another evidence that the 
Cooper pairs tend to gather between the ladders at small 
doping. 


V. DISCUSSION 

Let us first point out an important difference between 
the present results, obtained for weakly coupled ladders, 
and superconductivity in the hole-doped, square lattice 
Hubbard model. In the latter [18, 24] the order pa¬ 
rameter scales like J ~ At'^/U at large U. Here, A^, 
is nearly t7-independent in the range studied, as shown 
on Fig. 9. Changing the ratio U/t is typically accom¬ 
plished by applying pressure on the sample. However, 
in the case of Sri 4 _ 2 ,Caa;Cu 2404 i, changing the pressure 
would not only affect the value of t, but also of t' and, 
more importantly, doping, as carriers migrate between 
the chains and the ladders. Thus mapping a change in U 
to an experimentally accessible control parameter is very 
difficult. 

Let us now discuss the origin of the chiral supercon¬ 
ductivity that we have obtained in VCA. It is known 
that in the repulsive, large-!/ Hubbard model, the lattice 
symmetry and connectivity play an essential role in de¬ 
termining the symmetry of the order parameter. On the 
square lattice, real d-wave [dx' 2 _yi) symmetry fits well 
with the four-fold coordination. On the other hand, for 
triangular and honeycomb lattices, d + id ov chiral super¬ 
conductivity fits well with the three- and six-fold coordi¬ 
nation: Chiral states carrying a = 2 angular momen¬ 
tum avoid nodes in the order parameter bond values in 
real space, thus gaining condensation energy. Likewise, 
in an isolated and isotropic ladder, d-wave symmetry fits 
well because of the plaquettes. On the trellis lattice, we 
have elementary triangles, squares, and five-fold coordi¬ 
nation. A five-fold, odd number coordination in general 
accommodates a complex combination of d and s compo¬ 
nents. Clearly the amplitude of the chiral component of 




















superconductivity should increase with frustration, i.e., 
with t'. 

How to explain, then, that VGA and CDMFT disagree 
on the chiral nature of superconductivity? It may be that 
the small energy difference between the complex and real 
solutions shown on Fig. 6 cannot be resolved efhciently 
by CDMFT, but would be resolved if the same cluster- 
bath systems used in CDMFT were treated by Potthoff’s 
self-energy functional approach (a method called CDIA); 
in practice, this is impossible to do because of the large 
number of variational parameters involved. 

Despite this difference, the two approaches agree on 
important features: (1) The absence of superconductiv¬ 
ity at half-filling: the system is then a Mott insulator; (2) 
the existence of a ‘dome’ of d-wave superconductivity up 


to 15% to 20% doping; (3) the tendency of Cooper pairs 
to migrate from the inter-ladder regions to the plaquettes 
as doping is increased. A careful study of the order pa¬ 
rameter symmetry of the superconducting cuprate spin 
ladder compounds becomes important in the context of 
the possibility of chiral superconductivity found in this 
work. 
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